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Abstract 

We study the Cauchy problem with periodic initial data for the forward-backward heat equation 
defined by the J-self-adjoint linear operator L depending on a small parameter. The problem has 
been originated from the lubrication approximation of a viscous fluid film on the inner surface of the 
rotating cylinder. For a certain range of the parameter we rigorously prove the conjecture, based on 
the numerical evidence, that the set of eigenvectors of the operator L does not form a Riesz basis 
in C 2 {— 7r,7r). Our method can be applied to a wide range of the evolutional problems given by 
PT— symmetric operators. 



1 Introduction 

Analysis of the dynamic of a thin film of liquid which is entrained on the inside of a rotating cylinder 
is of great importance in lots of applications. For example when liquid thermosetting plastic is placed 
inside a rotating mould the best quality can be achieved if distribution of the liquid is as uniform as 
possible. More details about this application can be found in [26]. The same problem arises in the 
coating of fluorescent light bulbs when suspension consisting of a coating solute and a solvent is placed 
inside a spinning glass tube. The model for the coating was described for example in [7]. 

The lubrication approximation is used extensively to study flows in thin films. Under the assumption 
that the film is thin enough for viscous entrainment to compete with gravity, the time evolution of a 
thin film of liquid on the inner surface of a rotating in a gravity field cylinder can be described by the 
forward-backward heat equation: 

h t + Lh = 0, 0e(-7r,7r), te(0,T), (1.1) 

where 



Lh = ed e (sm 6 h e ) + h e , h{— k) = h(ir), e > 0. (1.2) 
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The effect of the surface tension is neglected in this linearized model derived by Benilov, O'Brien and 
Sazonov in [2]. 

We prove that the related to this equation Cauchy problem 

h\ t =o = h (9), h{-ir, t) = h(ir, t). (1.3) 

does not have a weak in the Sobolev sense solution h(6,t) even locally in time if ho(8) belongs to the 
class of finitely smooth functions with supp/io D (5, it — 5) ^ 0. 

The statement above can be roughly understood from the classic theory of the parabolic equations 
that states that regularity of a generalized solution depends on the regularity of the equation coefficients 
( in our case all coefficients are in C°°(— it, it)) and from the time-reversibility of the equation, i.e 
simultaneous change of the time variable t to —t and the space variable x to — x leads to the same 
partial differential equation. Time-reversibility and infinite regularity generally imply ill-posedness. 

The physical explanation of this explosive blow up of solutions is related to a drop of fluid that can 
be detached from the film in the upper part of cylinder, where the effect of the gravity is the strongest 

m p- 2i7]. 

The eigenvalues of the operator L were studied asymptotically, with application of the modified WKB 
approximation and numerically, with application of the analytic continuation method, by Benilov, 
O'Brien and Sazonov [2] and they came to a very interesting set of hypotheses: all eigenvalues of 
the operator L are located on the imaginary axis, they are all simple and the set of eigenfunctions is 
complete in £ 2 (—ir,ir) that is not typical for the ill-posed time-evolution problem. 

The analysis of the spectral properties of this operator was continued by Chugunova, Pelinovsky 
[9] and by Davies [TT]. Using different approaches they analytically justified that if the parameter 
|e| < 2 then the operator is well defined in the sense that it admits closure in £ 2 (—tt,it) with non- 
empty resolvent set. Analyzing tridiagonal matrix representation of the operator L with respect to 
the Fourier basis, Davies [TT] showed that L admits an orthogonal decomposition with respect to three 
invariant subspaces "H 2 '°(D), H 2 '°(C \ D), jmd Ker(L) = {c ]\, c G C} (see Section [4] below) and used 
this fact to prove that the nontrivial part L := L \ 7^ 2,0 (B) 7^ 2 '°(C\D) of L has a compact inverse of 
the Hilbert-Schmidt type. Therefore the spectrum of the original operator L is discrete with the only 
possible accumulation point at infinity. 

Under the additional condition (1/e is not integer) it was proved by Weir [30J that if there exists 
an eigenvalue A of the operator L, then /i = i 2 ^- is an eigenvalue of some symmetric operator, hence A 
can be only pure imaginary. The elegant proof is based on the continuation of the eigenfunctions into 
the Hardy space Ti 2 (H>) in the unit disk H> = {z £ C : \z\ < 1}. In this paper we sharpen this result 
showing that the additional condition above can be omitted. Note also that the similar problem was 
studied in the recent preprint [5] for a class operators that includes the operator L. 

It was shown numerically [91 [11] that the angle between the subspace spanned by the iV-first eigen- 
functions and the (N + l)-th eigenfunction of the operator L tends to as TV goes to infinity. This 
gives the numerical evidence that the eigenfunctions do not form a Riesz basis in £ 2 {—tt,tt) because 
the related projectors are not uniformly bounded. One of the main goals of this paper is to prove 
analytically this numerical conjecture justifying that the operator L is not similar to a self-adjoint. 

We also prove that the system of eigenvectors L is complete in £ 2 {— vr, it) (see Theorem 14.31 and 
Remark 15. 2D . Hence, this implies that L has infinite number of pure imaginary eigenvalues that 
accumulate to ±zoo. As a consequence, due to the linearity, the original Cauchy problem has infinitely 
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many global in time solutions which are linear combinations of harmonics e lXnt u\ n (x) where i\ n is an 
eigenvalue of the operator L and u\ n (x) is the related eigenfunction. 

The operator L is J-self-adjoint in the Krein space with J(f(9)) = f(ir—0) and therefore it belongs to 
the class of PT— symmetric operators. Interesting development of the spectral theory of PT— symmetric 
operators which are not similar to self-adjoint ones can be found in [20\ [2^1 [25j [5]. 

Notations: In the sequel, Ci, C2, ■ ■ ■ denote constants that may change from line to line but remain 
independent of the appropriate quantities. We also use h! , dgh, and hg for The symbol \\ denotes 
the function that identically equals 1 for 9 G (— tt, tt). Let T be a linear operator in a Hilbert space 
H. The following classic notations are used: Dom(T), Ker(T), Ran(T) are the domain, the kernel, 
and the range of T, respectively; cr(T) and p(T) denote the spectrum and the resolvent set of T; cr p (T) 
stands for the set of eigenvalues of T. We write f(x) X g(x) (x — > xo) if both f jg and g/ f are bounded 
functions in a certain neighborhood of xq. By D = {z G C : \z\ < 1} we denote the open unit disc in C. 



2 Analysis of the differential equation. 



One of the most general linear second-order differential equations with periodic coefficients that can 
be solved by a trigonometric series with a three-term recursion relation between the coefficients was 
studied by Magnus and Winkler p] and has the form 

(A + £cos(20))^f + Csm(20)^- + (D + Ecos{29))y = 0, (2.1) 
a9 z ad 

where A, B, C, D, E are constants. Under the additional condition that the coefficient A + B cos(29) 
does not have zeros located on the real axis they studied existence of the periodic solutions to (|2.1|) . 

In this section we study the basic properties of the differential equation £[h](x) = f(x) where the 
differential expression h is given by 

^=4H<)+^ 6 e <-*•*>• <2 - 2 » 

and use these properties to define the maximal periodic differential operator associated with £ and its 
inverse at the end of the section. This equation can be transformed to the form (|2.ip but all singularities 
are located on the real axis so the additional condition of Magnus and Winkler is not satisfied. 

Let / G C 2 {— vr,7r) and e > 0. Denote 2+ = (0, tt), 2- = (— tt, 0). Consider the differential equation 

£[h]( x )=f(x) a.e. on (-vr,7r) (2.3) 

assuming that 

the functions h and esin(0)/i + h are in AC\ oc {I- U (2-4) 
i.e., are absolutely continuous on each closed subinterval of X_ UZ+. 

Lemma 2.1. Let h satisfy (|2.4|) . Then h is a solution of the equation £\h\(x) = f(x) if and only if h 
has the form 



h{9) 



|cot(0/2)| 1/£ f kf - ( /(t)|tan(t/2)| 1/e dt J + I f(t)dt + kf, 9 G T±, (2.5) 

V J±ir/2 J JO 



where kf and k^ are arbitrary constants. 
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The proof is based on direct calculations. 



Proposition 2.1. Assume that e G (0,2). A function h G C 2 {— vr,7r) satisfies (|2.4p and is a solution 
of the equation £[h] = f with f G £ 2 (— 7r, 7r) i/ and on/?/ has the form 



h(6) 



cot — 
2 



l/e W9 



/(*) 



t 

tan — 
2 



dt+ f(t)dt + kf, 6 el±, 
Jo 



(2.6) 



where k 1 are arbitrary constants. 

Proof. Assume that h is an C 2 (— ir, 7r)-solution of £[h] = f. Then it has the form (|2.5p . Note that 



/(*) 



t 

tan — 
2 



l/e 



for any 5 G (0,tt). 



Therefore there exist the finite limit 





tan — 


l/e \ 




dt 


V J±tt/2 


2 





If Ci ^ 0, then ([23]) implies that |/i((9)j > |C 2 | fl" 1 / 6 for > small enough, where C 2 > 0. Since 
e < 2, we see that /i G" £ 2 (0,7r). This shows that C\ = 0, and therefore h has the form (|2.6|) on (0, 7r). 
Similarly, one can show that /i has the form (| 2 . 6 1) on (— vr,0). 

Let us prove that any function h of the form (|2.6p belongs to £ 2 (0, 7r) (the proof for C 2 (— ir, 0) is the 
same). It is enough to check that h G £ 2 (0, S) and h £ £ 2 (ir — 5, ir) for sufficiently small J > 0. 

For 6 G (0, o") , we have 



l/(*)l 



t 

tan- 



* < / t 2/e di =^3||/||£^ 1/2+1/£ . 



Hence, 



cot 



l/e 



tan ■ 



l/e 



dt < 2~ 1/e Cs 



(2.7) 



and we finally see that h G £ 2 (0, 5) 
For G (7T — 5, 7r), we have 



Hence, 



tan — 
2 



di<C 4 + 2 1 / £ ||/|| £2 H (7r-t)- 2 / £ dtJ <C 5 + 2 1 / 



£ 2 l 71 " 



vl/2-l/e 



cot • 



l/e r e 



\f(t)\ 



tan ■ 



l/e 



dt < c 6 {tt - ef/ s + evil/Ik^ 



a/2 



So /i G £ 2 (tt-o». 



□ 
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In particular, we have proved that 

lim h{9) = kf and lim h{9) = kf + f(t)dt 

8^±0 0— >±7t=f0 Jq 

hold for any £ 2 -solution h. This implies that the condition 

h is continuous on [— w, tt] and periodic (2-8) 

is fulfilled exactly when 

k+ = fcf and / ±]\ . 

Let the symbol W2 P (— vr,7r) stand for the subspace of the space W^— vr,7r) consisting of periodic 
functions, i.e., functions satisfying the conditions u^\tt) = u^(-tt) (i = 0, 1, . . . , k — 1). The norm in 
this space coincides with that of the Sobolev space W^(— vr,7r). 

Proposition 2.2. Let e G (0,2), / G £ 2 (-vr,7r), and J^f(9)d9 = 0. Then an C 2 -solution of £[h] = f 
satisfies (12, 8h if and only if 



HO) 



cot ■ 



l/e r e 



f(t) 



tan • 



l/e 



dt+ f f(t)dt + k% for a. a. G (— k,tv) 
Jo 



(2.9) 



where k\ is an arbitrary constant. Moreover, any function h of the form (|2.9|) possesses the following 
properties: 

(i) h G AC[-tt,7t], h' G C 2 {— w,tt), and 

\\tiU2 <ir(|^| + ||/|| £ 2), (2.io) 

where K is a constant independent of f . 

(ii) sm{9)ti G AC[-tt,tt] and (sm(6)ti)' G £ 2 (-vr,vr). 

Proof. To show that /i' G £ 2 (— 7r, 7r), it is enough to prove that h! G £ 2 (0, 5) for any 5 > small enough. 
Since 



•— ?— cot 1 /- ^ / /(t)tan 1 /- G (0,5), 



it is sufficient to show that 



J S 9- 2 - 2 / £ (J* f{t)t l / £ dt^j d9 < d jf* \f{9)\ 2 d0, 



(2.11) 



for any / G £ 2 (0,<5). 

Denote := /(t)i 1/e . Then (|2"TT|) takes the form 

2 



(2.12) 
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This is a weighted norm inequality for the Hardy operator. Applying [2T] (see also [23] and references 
therein), we see that 



sup ( [ S r 2 - 2 / s dt] 1/2 ( f (e- 2 ^) 

Be[o,S) \Je J \Jo v ' 



1-2 



1/2 



< OO, 



(2.13) 



and therefore (12. 12H holds true. It is easy to see that the latter implies (|2.10p and statement (i) of the 
theorem. 

If / -I'll, then sn\(Q)h' + h G W\ p (—ir,ir) and, by statement (i), so is sin(9)h'. □ 

Introduce the space X 2 = {h G W^— 7r, 7r) : (sin 6) he G W^(— 7r,7r)} endowed with the norm 

Nil = \\ h e\\h(-n,n) + IK sin #)Mwi(-^)- 

Denote by X® the subspace of X2 comprising the functions h with the property J h(6) dO = 0. As 
a consequence of the definitions, we obtain that if h G X2, then the function (sinO)hg is absolutely 
continuous (may be after a change on a set of zero measure) and 

(sin0) h e \ e=0 = (sin 6) he\e=n = (sin 6) /ie|e=-7r = 0. 

Let the symbol C 2 (—tt,tt) stand for the subspace of C 2 (— n, it) comprising the functions / with the 
property £ v f(0)dB = Q. 

We write below a set of corollaries of Proposition 12.21 We also give the alternative prove of this 
result using the Galerkin method in Appendix [Al 

Denote by L the operator acting in C 2 (—tt,7t) and defined by 



Lf = £[f) for / G Dom(L) := X 2 . 

Clearly, Ker(L) = {c]\, c G C}. 

Let us denote by L the restriction of L on C 2 (— it, it) (= (Ker(L)) ) , 

L:=L\ £j(-7r, vr), Dom(L) := Dom(L) n C 2 p (-ir, vr). 



(2.14) 



(2.15) 



It follows from the remark after (|2.8p that RanL C C 2 (— it, it). So L is an operator in the Hilbert space 



c 2 J-^). 



To find the inverse operator L , let us symmetrize (|2.9p as 



h(0) 



cot • 



l/e r e 



fit) 



tan ■ 



l/e 



tan • 



1/^ 



dt + ki for a. a. G (— tt.tt). 



(2.16) 



Solving the equation (h, 1) = 0, we get 



2tt 



cot ■ 



l/e r 8 



fit) 



tan • 



l/e 



tan ■ 



(2.17) 



So for / G Ran(L), we have L" 1 / = /i with /» defined by (l2TT6D - (l2TT7j) . 
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Corollary 2.1. 

(i) The operator L defined by {2.11$ is a closed operator in C 2 (—tt,tt) (with the dense domain X 2 ). 

(ii) Its kernel Ker(L) is the one- dimensional subspace of constants {c]\, c € C}. 

(Hi) The range Ran(L) of L is the orthogonal complement to Ker(L), Ran(L) = C 2 (— n, tt). 

(iv) The operator L defined by 112.15]) (and acting in Cp(—-K,ir)) has a compact inverse L^ 1 . 

(v) The operator L : X® — > Cp(—TT,ir) is an isomorphism of X® onto C 2 (— w, tt). 

Proof, (i)-(iii) Let k\ = k±(f) be the linear functional defined by (|2.17p . It is easy to see that 

k\ is bounded on C?(— tt,tt). (2-18) 

This and Proposition 12.21 imply immediately that Ran(L) = C 2 (—tt^tt). So Dom(L) = C 2 (—tt,tt). It 
follows from (|2.10j) and 

ll^'ll£2(-7r,7r) ^ II^II^^tt), h G £ p (-7T,7r), 

that L~ 1 is a bounded linear operator on Cp(— vr,7r). Therefore L is closed and so is L. 

(iv)-(v) It follows from (pHO]) and (pT[8j) that L _1 is bounded as an operator from £p(—7r, tt) onto 
X®- This proves statements (iv) and (v). □ 

Remark 2.1. The adjoint operation £* is given by 

It is easy to see that the adjoint operator L* is unitary equivalent to L and has the form 

L*h = f[h], Dom(L*) = X 2 . 

Indeed, L = JL*J, where (Jh)(6) = h(n — 9). In turn, this implies that the statements analogous to 
that of Corollary \2.1\ are valid for L* . Note that J is a signature operator, i.e., J = J* = J~ x , so the 
operator L is J -self-adjoint and belongs to the class of the PT -symmetric operations. 

3 The ill-posedness of the Cauchy problem for the forward-backward 
heat equation. 

The linearized model of the thin film dynamic (II, lj) was derived without taking into account the 
smoothing effect of the surface tension. It's very natural to expect that a drop of fluid will detach itself 
from the "ceiling" of the rotating cylinder and it will inevitably fall down that perfectly fits into the 
ill-posed nature of the Cauchy problem for the forward-backward heat equation (jl.ip . The intuition 
based on the classic theory of backward heat equation, says that global in time classic solutions can 
exist only for some class of analytic in vertical strip (— tt, tt) functions with exponentially fast decaying 
Fourier coefficients. 
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From there on in this section we assume that the parameter < e < 2. Let us consider the parabolic 
problem 

ht + Lh = 0, 0€(-7r,7r), t€(0,T), (3.1) 

h\ t=0 = h {9), u(-7r,t)=u(n,t). (3.2) 

Note that, after the change of variables 9^—9, equation (|3.ip can be replaced with the equation 

ut — Lu = 0. 

Let Q = (— tt, tt) x (0,T). We prove in this section that the problem (|3.1| . (|3.2p is ill-posed in the 
classes of finite smoothness. In what follows, the symbol (•, •) stands for the inner product in the space 
C 2 in the corresponding domain (Q, Qi, . . .). 

Definition 3.1. By a generalized solution to the problem \3. A3.£\) from the space ^2p(Q) we imply 
a function h G £ 2 (0, T; W^ p (-7r, vr)) n C([0, T]; £ 2 (-vr, vr)) such that 

h (e)v(e,o)de 

-TT 

for alive £ 2 (0,T;W^(-vr,7r)) : v t G C 2 (Q),v(9,T) = 0. 

Note that we can use other definitions of generalized solutions (see |19j). In particular, it is possible 
to prove that the above definition is equivalent to the following definitions. 

Definition 3.2. A function h G £ 2 (0, T; W^ p (-7r, tt)) n C([0, T]; C 2 {-tt, tt)) such that h\ t=0 = u (8) is 
called a generalized solution to the problem 113.1]) . jff.jj)) if 

-(h,v t ) - (sin 8 h ,vg) + (hg,v) = 

for all v G £ 2 (0, T; W^(-vr, tt)) : v t G C 2 (Q), v(0, T) = 0, v{9, 0) = 0. 

Definition 3.3. A function h G £ 2 (0, T; W^ p (-vr, vr)) n C([0, T}; C 2 (-tt, tt)) such that h\ t=0 = h (9) 
and ht G £ 2 (0, T; WiL (— tt, tt)) is called a generalized solution to the problem \3. 1\) . \3. 2fy if 

ht + Lh = 

and this equality holds in the space C 2 (0, T; W^V L (— tt, tt)) (W 2p 1 (— tt, tt) is the negative space constructed 
on the pair WiL(— tt, it), C 2 (—tt,tt)). 

Theorem 3.1 (nonexistence). Assume that there exists k G N such that ho G W2„(— 7r,7r) ; supp/io G 
(0, 7r) ; and ho ^ W^"^— tt, tt). Then there is no a generalized solution from the space W^piQ) for 
problem HTTp . 131$ . 

Proof. Assume the contrary, i.e., that such a solution exists. Let supp/io C (5,tt — 5). Put Qo = 
(6,tt-5) x (0,T/2). We have 

/*7T 

-(h,vt)-(sm9h e ,v e ) + (h e ,v)= / h (9)v(9,0)d9 (3.3) 
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for all functions v from Definition 13,11 such that suppi; C Qq U {(#,0) : 9 G (8, ir — 8)}. Make the 
change of variables r = T/2 — t. Then the function h is a generalized solution of the equation 

h T - Lh = 

in Qq and h satisfies the condition 

h\ T=T /2 = h (0). (3.4) 

Here, we understand a generalized solution in the sense that the function h G £ 2 (0, T/2; Wjj, (8, it — 8)) Pi 
C([0,T/2];£ 2 (5,7r -<5)) is such that 

/"7T — 5 

- (h, v T ) + (sin 9 h 9 ,v e ) -(h e ,v) = - / h (9)v(9, T/2) d0 (3.5) 



for all v G £ 2 (0,T/2;Wl(5,ir - 8)) : v t G £ 2 (Q ), v(9, 0) = 0, u| fl=<5 = v\ g=n - S = 0. Now construct a 
function /ii being a solution to the problem 

h T - Lh = 0, /i| e=5 = h\ e=7T _ s = 0, /i t= t/2 = h (9), 

in the domain Qi = (8, n — 8) X (T/2, T). Since > 1, we see that a solution of this problem exists and 
belongs 
identity 



belongs at least to the space W 2 ,' 1 ^!) (see [191 Theorem 3.6.1]). This solution satisfies the integral 



f'7T — S 

-(h 1 ,v T ) + (fanOh ig ,vo)-(h ig ,v)= J h (9)v(9,T/2) d9 (3.6) 

for all v G £ 2 (T/2, T; Wl(8, ir - 8)) : v T G C 2 (Q 1 ),v(9,T) = 0, v\ e=5 = v\ e=7T -6 = 0. Put Q 2 = 
(8, 7T — 5) x (0, T). Summing (|3.5|) and (|3.6|) we conclude that 

- (h>2,v T ) + (sm9h 2 e,v 9 ) - (h 2 e,v) = (3.7) 

for all v G £ 2 (0, T; Wr,(8, ir — 8)) : v T G £ 2 (Q 2 ), u(5, T) = 0, v\ e=5 = v\ e=7T ^ 5 = 0. Here the function 
/i2 coincides with h for t G (0,T/2) and with h\ for i G (T/2,T). So the function /i 2 is a generalized 
solution of the equation 

h T - Lh = 

in Q2 in the sense of the integral identity (|3.7p . By |19[ Theorem III. 12.1], we have u G C°°((52)- 
Therefore, h 2 {9,T/2) = hi(6,T/2) = h (9) G C°°(<5,7r - 5). This contradicts to the fact that h £ 
W k 2 ;\-8,8). □ 

Remark 3.1. It is easily seen from the proof that it does not matter where the support of ho lies. The 
main condition is that ho ^ W2 +1 (<5, 7r — 8) for some 8 > 0. So the theorem can be strengthened. 

Theorem 3.2 (instability). Assume that the problem §3. ( 13. £?]) is densely solvable in the following 
sense: for ho from some set K of smooth functions that is dense in the space W| p (— vr,7r) (k > 2), 
problem §3.1\) . \3.2^ has a generalized solution h in the sense of Definition 13.11 Then there is no 
constant c > such that, for every generalized solution of the problem \3.1\) . \3. 2\) with initial value 
ho G K , the estimate 

WH&iQ) < c IIM W f p (-7T,7r) ( 3 - 8 ) 

holds. 
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Proof. We use the arguments of Theorem 13.11 Let us find a function ho G W| p (— it, it) such that 
supp/io G (0,7r), and /io ^ W^j"^— vr, 7r). Find 5 > such that suppuo C (5, vr — 5). Put Qq = 
(5/2,7r-5/2)x(0,T/2). ' 

Assume the contrary, i.e., that (|3.8j) holds. Construct a sequence /io n € K : \\ho n — ^o|lw. fc (-vr w) ~^ ® 
as n — > co. In this case, 

IMIw^M/^-a^) 00 as (3-9) 

Denote by /i n the corresponding generalized solutions to our parabolic problem. As in the proof of 
Theorem 13. 11 the change of variables r = T/2 — t, shows that the functions h n = h n (9, T/2 — t) satisfy 
the integral identity 

/•7T-5/2 

-(K,v T ) + (smeh ne ,v )-(h n e,v) = - / h 0n (6)v(6,T/2) dO (3.10) 

JS/2 

for all v G £ 2 (0, T/2; W*(5/2, vr - 5/2)) : «t G £ 2 (Q ), KM) = 0, V | e=(5/2 = ^ =7r _ a/2 = 0. 

Thus, the functions h n are generalized solutions to a parabolic equation in Qq. Using Theorems 
III. 8.1 and III. 12.1 and Theorem IV. 10.1 in [19], we obtain that the function h n (8,T/2 — r) is infinitely 
differentiable in Qq and the norm of this function in any Hoder space # 2+a ' 1+Q / 2 (Qi) with Qi = 
(5',7r — 5') x (e ,T/2) (5' > 0, Eq > 0) is estimated by some constant depending on 5',£o,a, and 
the norm ||/in||£ 2 (Q )- In view of 13.81 with h n ,ho n substituted for h,ho and the fact that the norms 
||^0n|lw fc (-7r ?r) are bounded, we can assume that this constant is independent of n. As a consequence, 
we have the estimate 

ll / in(6 , ,0)|| w fc + i (3(5/4j7r _3 (5/4) < c\\h n \\ C 2 {Qo) (3.11) 
where the constant c is independent of n. Comparing 13.111 13.91 we arrive at a contradiction. □ 



4 The completeness property for the operator L 

As in the previous section we restrict the parameter to the interval < e < 2. In this section, we 
prove that the system of all eigenvectors and generalized eigenvectors of the operator L is complete in 
C?{— 7r,7r). In particular, this implies that L has infinite number of eigenvalues. 

Denote by H 2 '°{B) and H 2 '°(C \ D) the subspaces of the Hardy spaces H 2 (B) and H 2 (C \ B) , 

respectively (see e.g. [121 Section 2.1]) that are orthogonal to the function "|| (= 1). In the sequel, we 

use the standard identification of the function u(z) G H 2 (IS) with the function u(e ) := lim u(re ie ), 

r—fl—0 

which belongs to C?{— 7T, tt), and also use the similar agreement for u(z) G H 2 (C \ D). Then Tt 2 (H>) 
and n 2 {C\B) are the subspaces of £ 2 (-tt, vr) and H 2 (U>) nH 2 (C\ D) ={c% c G C}. In these terms, 
the space £ 2 (— vr,vr) admits the orthogonal decomposition 

£ 2 (-tt, vr) = H 2 '°(B) e {c 1|, c G C} H 2 >°(C \ D). (4.1) 

Define the operator Lg n in the Hilbert space C 2 {— vr, vr) by L^ n h := £[h], Dom(Lsn) = Pfi n , where I 
is the differential expression defined in Section [2] and Pfi n is the set of finite trigonometric polynomials 

N 

h(9) = (2VT)- 1 / 2 v ne ind , N < oo, v n G C. 

n=-7V 



10 



It is easy to see that Lg n is densely defined, and hence the closure 

exists as an operator in C 2 (—tt,tt). 

Let L be the indefinite convection-diffusion operator defined by (|2. 14[) and let L be its restriction 
defined by (|2.15p . It is easy to see that Remark 12.11 implies that Lg n = L* and therefore L m \ a = L. 
Below we give another proof of this fact using the results of |11| . This proof allows us to use the 
orthogonal decomposition of L obtained in [IT] (see also [9]). 

Proposition 4.1 (Theorems 11 and 13 in [11]). 

(i) The operator -£/ mm admits the orthogonal decomposition L m { n = L_ © © L + with respect to fl^. 

(ii) The operators L± are invertible, and their inverses L^} are Hilbert- Schmidt operators. 
(Hi) L+ 1 and (— LJ)~ l are unitary equivalent. 

Proposition 4.2. If e <2, then L = L m i n and L = L_ © L + . 

Proof. By Corollary 12.1 1 (i). L is a closed extension of Lfi n . Hence L m \ n C L. Let us show that L = L mm . 

By Proposition 14.11 the operator L^ 1 © LZ 1 is compact and is acting in C 2 ( — it, it). So er(L_ © L + ) 
is at most countable. Then 

i?an(L_ © L + - XI) = C 2 p ( — 7r,7r) for any A £ p{L- © L + ). (4-2) 

By Corollary 12.11 (iv), the operator L" 1 acting in C? p {— vr, it) is compact. So L possesses the same 
properties, that is, o~(L) is at most countable and (|4.2p holds for L. 

Assume that L + © L_ £ L (which is equivalent to L m \ a S L). Then Eqs. (|4.2p for L + © L_ and L 
imply that A G o~ p (L) whenever A £ p{L) fl p{L + © a contradiction. □ 

Definition 4.1 (e.g. [H]). By & p , < p < oo, we denote the class of all bounded linear operators A 
acting on a Hilbert space H for which 

\ A \p ■■= (^t sP M^J <0 ° 

where Sj(A) are singular numbers of A, i.e, eigenvalues of the self-adjoint operator {A* A) 1 / 2 that are 
enumerated in decreasing order, counted with multiplicities. 

Two S-classes were given special names: 62 is the class of Hilbert- Schmidt operators and Si is the 
class of nuclear operators. It was proved by Davies [11] that the operators A± , where 

A ± := -iL±, (4.3) 

belong to the class ©2 and so does L _1 . We prove in this section that, actually, the operator is 
nuclear. 

We need the following result (see [13], its weaker version can be found e.g. in [HI Section III.7.8]): 
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Theorem 4.1 (Gohberg, Markus). If0<p<2, then the linear operator A acting in a Hilbert space 
H belongs to 6 P if and only if for at least one orthonormal basis {ej} of H the inequality 

oo 

^L4 ej f<oo (4.4) 

3=1 

holds. In addition, 

oo oo 

3=1 j,k=l 



It follows from Theorem 11 and Eq. (15)] that the operator iL + (= A + ) in the Fourier basis 
{ e n}i°' e n(&) = e ind , is represented by matrix (p m ,n) which has the following properties 

\Pm,n\ < Ci m~ l+1 l e n" 1_1 / £ , m<n, 

\pm,n\ < Cim- 1 - 1 ' 6 n- 1+1 / e , n < m. [ ' ' 

Proposition 4.3. The operator L^ 1 belongs to the class &\ for any e G (0,2). More precisely: 

(i) if £ G (0, 1], then L^ 1 G 6 p for any p > 2/3, 

(ii) if e G (1, 2), then Z" 1 G & p for any p > 2e/(e + 2). 



Proof. 

(i) ££(0,1]. 



Using (|4.5p . one can obtain that 

oo n oo 

m=l m=l m=n+l 

/>oo 

< Cfn" 3 + C 2 n" 2+2/e / m- 2 - 2/e dm < C 2 n" 3 . (4.6) 

J n 

Hence, 

oo oo 

J2W L +en\\ P c i<C 3 J2n- 3p / 2 . 

n=l n=l 

So the Gohberg-Markus criterion (Theorem 14. ip shows that L^ 1 G S p for any 2/3 < p < 2. Thus L^ 1 
belongs to S p for any p > 2/3 and so does L^ 1 due to Propositions 14.11 (ii) and 14.21 
(ii) £G (1,2). 

Using (|4.5p . one can obtain that 



ll^enll^ < C 2 n- 2 ~ 2 / £ ^ m~ 2+2 / £ + C 2 n- 2+2 / £ ^ m" 2 ^ < Cfn" 1 - 2 /* + C 2 n' 3 < C z n~ x ~ 2 l\ 

m=l m=n+l 

Therefore, L^ 1 G & p whenever p(e + 2)/(2e) > 1. □ 
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Although the weaker result that the operator L € &„ for p > 1 can be obtained directly from 
the factorization of the operator L found by Chugunova and Strauss in [lOj . the fact that the operator 
L _1 actually belongs to the class of nuclear operators Si is crucial, as you see below, for the proof of 
completeness. 

Following [144 Section IV. 4], we will call an operator T acting in a Hilbert space H dissipative if 

Im(Tf,f)>0 for all f £ Dom(T). (4.7) 
Proposition 4.4. The operators L + , (— L + ) _1 , — L_ and LZ 1 are dissipative. 



Proof. Using the tridiagonal matrix representations of A + {= —iL+) with respect to the Fourier basis 
{e in0 }f° (see 0Q3]), we get: 



^4+ — (fln.m)* 



1 


— £ 








e 


2 


-3s 








3c 


3 


-6e 








6s 


4 



a n ,n = n, a n _i jn = ^n(n - 1), a„ >n+ i = -|n(n + 1), n=l,2, .... 



(4- 



This representation implies that 



Im(L+h, h) = Re(A + h, h)>0 



(4.9) 



for all h G W 2 '°(B) n Pfi n . Since L = Lfi n , one gets (I4.9P for all h <G Dom(L + ), i.e., L + is dissipative. 
Substituting /i = (L+)- 1 / into (IQ]) . we see that so is (— L + ) . Proposition 14.11 (iii) completes the 
proof. □ 

Theorem 4.2 (Lidskii, see e.g. Theorem V.2.3 [14] ). If the dissipative operator A acting in a Hilbert 
space H belongs to the class &\, then its system of all eigenvectors and generalized eigenvectors is 
complete in H . 



Now the main result of this section can be obtained using Propositions 14. lj, 14.31 14.41 and Lidskii's 
theorem. 

Theorem 4.3. The operator L has infinitely many eigenvalues. The system of its eigenvectors and 
generalized eigenvectors is complete in C 2 (—7r,ir). 

Proof. By Propositions 14.31 14.41 and Lidskii's theorem, the system of all eigenvectors and generalized 
eigenvectors of the operators , LZ , and L _1 are complete in W 2 ' (D), 7Y 2 '°(C\D), and £?(— vr,7r), 
respectively. Since Ker^L" 1 ) = {0}, all generalized eigenspaces of L^ 1 corresponding to its eigenvalues 
a n are generalized eigenspaces of L corresponding to eigenvalues A n = l/a n . Since all eigenvalues of the 
compact operator L _1 have finite algebraic multiplicities, we see that L _1 and L have infinitely many 
eigenvalues. The finite-dimensional spectral mapping theorem implies that the completeness property 
holds for L and, consequently, for L. □ 
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5 Pure imaginary eigenvalues and the Riesz basis property 



In this section, the following result is obtained: all eigenvalues of L are pure imaginary (this statement 
was proved by Weir [30] under the additional assumption that 1/e ^ Z). We use this fact to prove that 
eigenvectors of L do not form a Riesz basis in C 2 (— vr,7r). 

Recall that L + is an operator in the Hilbert space H 2 '°(B) = H 2 (I}) G {c ]\, c £ C}. We identify the 
function u(z) G H 2 (D) with u(e w ) G £ 2 (-7r,vr). Note that 

u ±.]\ is equivalent to u(0) = 0. 

So the last equality holds for all u G H 2 '°(D). 

Let u(z), z G D, be an eigenfunction of the operator ^4+(= —iL + ). Consider the restriction u of the 
function u on the interval [0, 1) C D, 

u(x) = u{x) for i € [0, 1). (5.1) 

The following proposition obtained by Weir [30J shows that if u is an eigenfunction of the operator 
A + {= — iL + ), then its restriction u is a solution of a Sturm-Liouville eigenvalue problem with real-values 
coefficients. 

Proposition 5.1 ([30]). Assume that A is an eigenvalue of the operator L + , u(e ld ) is a corresponding 
eigenvector, and u is its restriction defined by ( f5.il) . Then b[u](x) = fj,u(x) for all x G (0, l)(c D), 
where fj, = —2i\/e and the differential expression b is defined by 

b[u] = --(puj, 
w 

P (X) = (1 - X) 1+1/E (X + l) 1 - 1 ^, W(X) = 2^(1 " X) 1/£ (X + 1)- 1/E . 

Let -B max be an operator in £ 2 ((0, 1); w) associated with the differential expression b[-] and defined 
on its maximal domain 

B max u=b[u], MB max ) = {uG£ 2 ((0,l);t«):u,u'GACi oc (0,l), b[u] G £ 2 ((0, 1); w)}. 

Note that all points of the interval (0, 1) are regular for the differential expression b, but the endpoints 
and 1 are singular (1 is singular since pT x £ £ 1 (l/2, 1)). 

Proposition 5.2. Let e > 0. 

(i) b is in the limit-point case at 0, 

(ii) b is in the limit-point case at 1 exactly when e < 1. 

(iii) -B max is self-adjoint in £ 2 ((0, 1); w) exactly when < e < 1. 

Proof (i) Clearly, "|| is a solution of b[u] = and "||^ £ 2 ((0, 1/2); Weyl's alternative (see e.g. [29l 
Theorem 5.6]) completes the proof. 
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(ii) The general solution of b[u] = on (0, 1) is 

f x 1 

u(x) = ki / — — ds + k2, fci,^eC. 
Ji/2P(s) 

Clearly, 

u(x) x ki / (1 - s)~ 1_1/£ (is + A: 2 x fci(l - x)~ 1/e + fc 2 , x 1 - 0. 

A/2 

Hence all solutions of b[u] = belong to £ 2 ((l/2, 1); w) if and only if e > 1. 

(iii) follows from (i) and (ii). □ 

Proposition 5.3. If u{e ld ) G W 2 '°(D) ; £/ien i/ie function u defined by 115. 1\) belongs to £ 2 ((0, l);w). 



Proof. By ([5]), we have u(0) = 0. Since u is analytic at 0, we see that u G £ 2 ((0, 1/2); w). The measure 
w(x)dx on [1/2, 1) induces a measure M(S) := / Sn [i ^ w{x)dx on D. For any sector 

5 = {re ld :l-l<r<l,\9-9 \<l}, I G (0, 1) 



we have 



Af (5) < 2 / dx < 21 

•/max-fl/2,1-0 



c{l/2,l-Z} 

since rnax a . e [ 1 / 2 ,i) to(a^) < 2. So M(-) is a Carleson measure (see e.g. [121 Sec. 4.3]). Therefore, 

u(x)\ 2 w(x)dx < Ci\\u\\^2o , 

1/2 

where C\ is a constant independent of u. This completes the proof. □ 

Proposition 5.4. Lei e G (0,2). Then all eigenvalues of the operator A+(= —iL+) are real and 
positive. 



Proof. Let u{e ie ) G H 2 '°(B) and L + u = An. By Proposition E3 u G £ 2 ((0, 1); w). Proposition O 
implies that b[u] = uu with /i = —2iX/e. Let us split the interval (0,2) into two parts. 

If e < 1 then the proof is simple. Clearly, u G Dom(B ma _ x ), and therefore u is an eigenvalue of the 
nonnegative self-adjoint operator B m&x . Thus, u > 0. 

If 1 < £ < 2 then the proof requires additional analysis. By Proposition 12.21 9( eid ) := ^f( eie ) ^ 
£ 2 (— vr,7r). It is easy to see from the representation u{e %e ) = (27r)~ 1//2 XmLi v n^ me and that g{e %e ) G 
W 2 (D) (on the other hand, the latter follows from Theorem 16]) and g(e ld ) = lim r ^i-o g(re ) 
where g(z) = z^jf-, z G D. By [12, Problem II.5 (a)], |g(x)| < IMI^l - |^| 2 ) _1/2 for x G (0,1) and 
therefore, for x G (1/2, 1), 



du(x) 



dx 



< Ibll^lxr^l - |x| 2 )~ 1/2 < Ci(l - x)- 1 / 2 . (5.2) 



By |29} Theorem 5.8 (ii)], the operator Bm denned by B-uu := b[u] on the domain 

Dom(Bi\) := {u G Dom(B max ) : []\,u]i = 0}, [1|,u]i:= limp(x)u'(x) 
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is self-adjoint. Note that the limit ["||,u]i exists for any u £ Dom(B max ) due to |29t Theorem 3.10]). 

It follows from (|5.2p that, for any eigenvector u(e l8 ) of A+, its restriction u belongs to Dom(Bn). 
Indeed, it was shown in the step (1) of the proof that u S Dom(B max ). On the other hand, it follows 
from that 

[1|,u]i = lim (l-x) 1+1/£ (x + l) 1 - 1/E -u'(x) = 0. 

x— »1— 

So /i is an eigenvalue of the operator B^ = BT,. 

It follows from (|5.2|) that u(x) = u(l/2) + J"^ u'(t)dt has a finite limit as x — > 1 — (this fact also 
follows from Theorem 16]). Therefore, 

«)jca((o,i)^o(*)) = -J {p{x)v!{x))'u(x)dx = J p(x)|u'(x)| 2 dx-^hm o p(x)u / (x)u(x) (5.3) 

= I p(x)|u'(x)| 2 dx > 0. (5.4) 
J o 

Thus, Bm > and therefore fi > 0. 

Finally, note that Ker{L + ) = and therefore (j, ^ 0. □ 

Remark 5.1. For e £ (0,2) suc/t i/iai 1/e Z, a slightly different form of Proposition \5.4\ was 
proved in jl30[ Theorem 2.3] by means of Proposition HOI the Frobenius theory, and the deep analysis 
of eigenvectors of the corresponding recursion relation given in Sections 2 and 3]. Our proof that 
removes the condition 1/e ^ Z, is based on the description of Dom(L) given in Section\M and the Hardy 
spaces theory. 

Theorem 5.1. If e € (0, 2), then the set of eigenvectors of the operator L does not form a Riesz basis 
in C?{— 7T, 7r). 

Proof. Assume that the set {u n }^° of all (linearly independent) eigenvectors of L form a Riesz basis in 
C 2 {— vr,7r). Then Proposition 15.41 implies that iL is similar to a certain self-adjoint operator Q. That 
is, there exists a bounded and boundedly invertible operator S such that S Dom(Q) = Dom(L) and 
iL = SQS- 1 . 

The spectral theorem for a self-adjoint operator implies that, for arbitrary uq £ Dom{L){= X2), the 
problem 

ut + Lu = 0, u \t=o= uq, iGl, 

has a unique solution u(-,t) in the sense of [16, Definition 1. 1.1 and Eq. (1.1.2)] (such solutions are 
sometimes called strong solutions). Moreover, this solution has the form u(-,t) = Se~ %t ® «S _1 iio(-)- 
Therefore, for any T > 0, 

u G C([0,T];X 2 ) c C([0,T];W^(-7r,7r)), 
ut(;t) € C([0,T];£ 2 (-7r,7r)), and Lu{-, t) € C([0, T}; £ 2 (-vr, tt)). 

It is easy to see that u is a generalized solution of (|3.ip , f|3.2|) in the sense of Definitions I3.HI3.31 Since 
e -itQ j g a un j^ ar y operator, 

ll U (-^)ll£ 2 (-7r,7r) < \\S\\ WS^W \\Uo\\c2(-it,k), * S M. 
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Hence, for any T > 0, we have 

f T 

J o IM> ^ CT W U 0\\%(-Tr,n)> 

where C = \\S\\ WS^ 1 1| 2 < oo. The latter contradicts Theorem [221 since Dom(L) is dense in W 2 ^— vr, n) 
and ||u || £ 2 ( _ 7ri7r) < IKollw^C— 7r,7r)- D 

Remark 5.2. Arguments of Proposition [5^] show that all eigenvalues of the operators L± are simple. 
Combining this with Theorems \4-3\ and \5.1\ one can show that (under the assumption e G (0, 2)) 
the system {u^f 3 of eigenvectors of L is complete in C 2 {— n,ir), but does not form a Riesz basis in 
£ 2 (-7r,vr). 

Indeed, let us show that all eigenvalues of L + are simple in the case e G (1, 2). Assume thatuxie 16 ) G 
H 2,0 (n) \ {0} is a generalized first order eigenvector of L + , i.e., (L + — XI)u\ = u, where u(e lS ) G 
7^ 2 '°(D) \ {0} and L + u = An. Consider the restrictions u, iq of the functions u and u\ on the interval 
[0,1) CD. It follows from the proof of Proposition \5J\ that u, Ui G Dom(B^) C £ 2 ((0, 1); ui). On the 
other hand, computations analogous to that of [30, Lemma 2.1 and Theorem 2.3] show that b[ui] — uui = 
— y-u with /i = — Therefore ui is a generalized eigenvector of the self-adjoint operator Bm, a 
contradiction. The proof for the case e G (0, 1] is similar. 

We would like to note that the linear partial differential equation is an interesting example 

when the nature of explosive blow-up and instability of solutions has its roots not in location of the 
eigenvalues but in geometric properties of the eigenfunctions. 



6 Further discussion 

When eigenfunctions related to neutrally stable eigenvalues of some linearized problem form the com- 
plete set, representation of a solution of the nonlinear problem as a series of these eigenfunctions is 
one of general approaches to the nonlinear stability problem. The lack of a basis property of the 
eigenfunction set is an obstacle for the applicability of this particular method. 

Due to the ill-posed nature of the forward-backward heat equation all eigenmodes are linearly un- 
stable [2] and it is common to use the smoothing effect of the surface tension to stabilize them . The 
lubrication approximation that takes into account the influence of the capillary effects and/or surface 
tension leads to the initial value problem for the fourth order nonlinear partial differential equation 
[8] . Some stability properties of its linearization were studied in O HI [6] . They came to the conclusion 
that almost all but some first modes are getting stable even if the surface tension is relatively weak. 

We would also like to mention that the main assumption about the parameter range |e| < 2 comes 
naturally from the theory of mixed type equations and for the case when |e| > 2 all properties of this 
backward forward heat equation can be changed significantly. 

Let us consider the equation 

k(x,t)u tt + a(t,x)u t + Au = 0, x G tt, t > (6.1) 

where the coefficient k(x, t) can change sign in the domain where the operator is considered. So equation 
(16. ip is an equation of the mixed type, i.e. it is of the same type as the well-known Tricomi equation. 



17 



On the lateral boundary of the cylinder fi x (0, T) we pose the Dirichlet boundary condition and there 
are two additional boundary condition on the lower and upper base of the cylinder: 

u\ t=T = 0, ut\ s+ =0, u t \ s -=0, S + = {(0,x) : k(x,0) > 0}, S~ = {(T,x) : k(x,T) < 0}. 

This boundary value problem and close problems were studied by many authors (see, for instance, 
|27t [2"B] ) . It was demonstrated that the condition 

a-^>5 >0 V(z,t), 

where 5q is a positive constant, ensures the existence of generalized solutions to the above-described 
boundary value problem. Stronger conditions of the type 

B _NPj -!)>«,>„ 

ensure existence of smooth solutions and uniqueness of generalized solutions. The existence of solutions 
of non-linear forward-backward heat equations was studied by Hollig [15J and by Pyatkov [22J. Among 
last results devoted to the nonlinear forward-backward parabolic problems we would like to mention 
Kuznecov papers jTTl [18] . 

Acknowledgement. The authors thank A. Burchard and E.B. Davies for useful comments and 
discussions. The research of M. Chugunova is supported by the NSERC Postdoctoral Fellowship. 
I.M. Karabash would like to thank P. Binding for the hospitality of the University of Calgary. 



A Proof of Corollary 12.11 using Galerkin method. 

The second proof of Corollary 1 2. 1\ Let {uij}JL l be a basis for the Hilbert space H = {h € W| p (— vr, n) : 
J-n h{0) d6 = 0}. Find functions ipj such that ipj$ = Uj, J_ d9 = 0. We look for an approximate 

solution to equation ()2.3f) in the form 

n 

h n — ^ ^ Cj n ifj, 
3=1 

where the constants Cj n are determined from the system of algebraic equations 

(Lh n ,Uj) = (f,uj), j = l,2,...,n, (A.l) 

(the brackets denote the inner product in £ 2 (— vr, it), i.e., (h,v) = h(9)v(8) d6). 

Multiply (jA.lh by c, n and summarize the equalities obtained. We arrive at the relation 

(Lh n , h n e) = (/, h n g). 

Integrating by parts we derive the estimate 

\\h n e\\c 2 (-TT,TT) < c ll/ll£ 2 (-7r,7r)) ( A - 2 ) 

where c is a constant independent of n. This estimate implies that the system (lA.ip is solvable. Note 
that there exists a constant c\ independent of n such that 

IIMl£ 2 (-7r,7r) ^ C \\hne\\c 2 (-n,n) ( A - 3 ) 
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From (|A.2p . ()A.3|) we conclude that there exists a subsequence h nk and a function h £ W 2 (— 7r,7r), 
h(—ir) = /i(vr) and f_ /i(#) d# = 0, such that 

h Uk ^h in £ 2 (-tt, ir), h nk e ^ he weakly in £ 2 (-7r, 7r). (A. 4) 

Multiply (lA.lj) with n = by constants a.,- (1 < j < m < n^) and sum the results. Fix m assuming 
that nfc > m. We infer 

m 

-e(sin 6h nk o,uj e ) + {h nk6 ,u)) = (f(9),u>), u = ^ccjUj. 

3=1 

Passing to the limit as k — > oo we arrive at the equality 

- e(sin6 , /i6»,W6i) + w) = (f(0),u), u = y^ajUj. (A. 5) 

i=i 

m 

The functions lo of the form u = ^ Oj^j are dense in and thus (|A.5[) holds for all functions in ii". 

i=l 

Due to our condition for the function /, we can see that (|A.5P also holds for all functions of the form 
uj + c (c is an arbitrary constant) and therefore for all functions in W|p( — 7r,7r). In particular, it holds 
for u; E 7r, 7r). From the definition of the generalized derivative (in the Sobolev sense) we have 

that there exist the generalized derivative (sin Q}iq)q and 

e(sm9h 9 )e = (/- M G £ 2 (-7r,7r). 

Thereby, sinOhg £ W 2 (— 7r,7r). Integrating by parts in (|A.5|) we obtain that the equation (|2.3p is 
satisfied almost everywhere on (— 7r, 7r). We have proven that L is an isomorphism of H® onto £^(— vr, 7r). 
The remaining assertions more or less obvious. 

Remark A.l. We can take the set {sin j9, cos jO}'jL 1 rather than an abstract basis {uJj}. 

In view of the embedding theorems, the operator L _1 : £ 2 (— 7r,7r) — > £ 2 (— vr,7r) is compact and thus 
the spectrum of L is discrete with the only accumulation point oo. Assume the contrary that there 
exist AsM such that 

Lh = Xh, he H%. 

Multiply the equation by h and integrate the result over (— 7r,7r). Integrating by parts and taking the 
real part, we arrive at the inequality 

II M 2 ^.^) < o. 

We have used the equality Re(h, ho) = 0. Thus, h = 0. □ 

Theorem A.l. Under the condition 1 — e(k + l/2) > 0, for every f G W 2p (— tt,7t) there exists a unique 
solution to the problem (2.1) such that u E W 2 j~ (— 7r,7r) ; sin {—TT,ir). 



Proof. We can use the same arguments as in the second proof of Corollary 12. 1[ It is not difficult to 
show that there exist set of constants Aj > such that in the equivalent inner product in the space 
W 2 fc (-7r,vr) 

k 
3=0 



(M)* = I>^ (J V 
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we have the inequality 

(Lh,he) k > 5 \\h\\ 2 w , +1{ _ n ^ y Vh G W 2 fe p +2 (-7r,vr), 

where the constant 5q > is independent of h. Next we apply the same arguments as those in the 
theorem 1 but we use the inner product (h,v)k rather than the inner product (h,v) in C?{— w, it). So 
the Galerkin method is applicable here. □ 

Remark A. 2. It is also possible to use some functional arguments based on the Hahn-Banach theorem. 
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